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The gauge-independent phenomenon of color confinement in Yang-MiUs theory manifests itself 
differently in different gauges. Therefore, the gauge dependence of quantities related to the infrared 
structure of the theory becomes important for understanding the confinement mechanism. Par- 
ticularly useful are classes of gauges that are controlled by a single gauge parameter. We present 
results on propagators and the color-Coulomb potential for the so-called A-gauges, which interpo- 
late between the (minimal) Landau gauge and the (minimal complete) Coulomb gauge. Results are 
reported for the SU(2) lattice gauge theory in three and four space-time dimensions. We investigate 
especially intermediate and low momenta. We find a continuous evolution of all quantities with the 
gauge parameter, except at zero four-momentum. 

PACS numbers: ll.15.Ha 12.38.Aw 14.70.Dj 



I. INTRODUCTION 

The investigation of the confinement mechanism in 
QCD and Yang-Mills theory has progressed during the 
last decade, but remains an unsolved and challenging 
problem (see e.g. [H for a brief review). A particularly 
important question is that of gauge dependence. 

The confinement phenomenon has been investigated in 
various gauges, with different degrees of sophistication, 
methods, and observables. Still, at the present stage it is 
not yet clear to what extent the confinement mechanism 
itself is gauge-dependent. It is, however, sure that its 
manifestation in the structure of gauge- dependent quan- 
tities, like correlation functions, is different in different 
gauges. 

For a full understanding of the confinement process, 
it is thus necessary to understand this mutability among 
different gauges. For such a task, classes of gauges are 
of special interest which offer parameters to deform the 
gauge condition smoothly. This allows one to track the 
changes of correlation functions along the gauge orbit. 

One such class are the A-gauges [2[ Q , defined by the 
gauge condition 



(1) 



where d is the number of space-time dimensions and the 
A^'s are a set of gauge parameters. E. g., the Landau- 
gauge case will be given by A^ = 1 for all fi. In princi- 
ple, it is necessary to add a prescription for dealing with 
Gribov-Singer effects [1, 0- We employ here the mini- 
mal extension of the perturbatively defined gauge ([1]) to 
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the non-perturbative domain by selecting a random Gri- 
bov copy among all Gribov copies belonging to the first 
Gribov region. This region is defined by those gauge 
copies for which the associated Faddeev-Popov operator, 
defined below in equation ([3]), is non- negative. 

Clearly, if one of the A^'s identically vanishes, the 
gauge condition is similar to a would-be Coulomb gauge, 
but has to be supplemented by further constraints to 
yield even a perturbatively complete gauge condition. 
In contrast, a complete Coulomb gauge condition is ob- 
tained by performing a limit A^ — >■ for a certain /j, Q . It 
is not clear, whether the different lattice definitions of a 
perturbatively complete Coulomb gauge 0] yield coin- 
ciding correlation functions. The results we present here 
indicate that this may be the case, at least for a subclass 
of these definitions. However, in case of sufficiently non- 
smooth completions of the Coulomb gauge some effects 
could be present 

Here, we will concentrate on Landau-like A-gauges, in 
which we choose only one of the A^'s different from 1. 
As the lattice calculations presented here will be per- 
formed in Euclidean space-time at zero temperature, the 
direction is irrelevant, but will conventionally be chosen 
to be the time direction. Therefore we set Ao = A and 
Ai = A2 = • • • = Xd-i = 1 henceforth. 

In these gauges the simplest correlation functions are 
the gluon propagator and the Faddeev-Popov-ghost prop- 
agator. We will use here lattice gauge theory to deter- 
mine these propagators in the case of SU(2) Yang-Mills 
theory. Furthermore, in Coulomb gauge there exists an 
additional color-Coulomb potential that bounds the con- 
ventional Wilson potential from above [1]. A confin- 
ing color-Coulomb potential is thus a necessary condi- 
tion for confinement, and we shall also investigate how 
this potential develops from Landau to Coulomb gauge. 
There exists an associated order parameter jthe so-called 
residual-gauge-symmetry order parameter , which sig- 
nals whether the residual gauge symmetry that exists in 
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Coulomb gauge is unbroken. Also its evolution with the 
gauge parameter A will be determined here. 

A word of caution is due here. The behavior of the cor- 
relation functions in the asymptotic infrared domain has 
not yet been completely explained even in Landau gauge . 
In fact, despite great effort invested in analytic |lC)l4l3l | 
and numerical |14l4l7| studies, a consistent description 
unifying the different approaches is still being elaborated. 
For recent reviews of analytic and lattice studies see re- 
spectively [l8l - l20j and Ho, [U. Since the gauge-fixing 
process in the interpolating gauges is even more compu- 
tationally expensive than in ordinary Coulomb or Lan- 
dau gauge, it will not be possible here to reach volumes 
permitting to answer the question whether a scaling-type 
or a decoupling-type behavior is observed for general A- 
gaugcs. The main interest is therefore on the general de- 
formation of the relevant quantities with the gauge pa- 
rameter. Note that functional-equation studies predict 
the existence of a scaling-type solution in these gauges 
in four dimensions [l^]. For three dimensions, this will 
be shown in the appendix. Following the lines of the ar- 
gument given in |18| . it is easily conceivable that also a 
decoupling type of solution exists in this class of gauges. 
In any case, we will not attempt to settle this question 
here. 

The lattice definition of the A-gauges and the details 
of the production of the corresponding configurations are 
outlined in Section|lIl The results for the propagators are 
presented in Section IIIII Definitions of the propagators 
are also introduced there. Results for three dimensions 
are given in Subsection IIII Al while those for four di- 
mensions are given in Subsection IIII Bl Definitions and 
results for the color-Coulomb potential and associated 
quantities will be discussed in Section HVl A summary is 
given m Section El 

The results presented in this paper extend earlier stud- 
ies [23} . In particular, we complement our previous in- 
vestigation [g , in which a qualitative difference of some 
parts of the gluon propagator was observed in d = 4. 
Some preliminary results were presented in Ref . (2^ . See 
also [26] for investigations of other setups and quantities 
in this type of gauges. 

II. DEFINITIONS AND SET-UP 
A. Generation of gauge-fixed configurations 

Minimal A-gauges on the lattice are a straightforward 
gen eralization of the conventional minimal Landau gauge 
[231 [23. In particular, to fix the gauge, a gauge trans- 
formation g(x) is determined, which minimizes the func- 
tional 

^^^^""^^ = 1 - ^ E E K(^) + • (2) 

Here a is the lattice spacing and V = iV^ the lattice 
volume, [/^t € SU{Nc), with Nc = 2 for the present 



SU(2) case, is a link variable in the lattice configuration 
of the gauge field, and its gauge transform Uj^ is given by 

Ujiix) = g{x)U^ix)g+ (x + ae^). 

For a more detailed discussion of this gauge condition, in 
particular for the limit of A — > and the connection to 
Coulomb gauge, see 

Minimizing the above functional guarantees that the 
gauge-fixed configuration is within the first Gribov hori- 
zon, defined by the first null eigenvalue of the lattice 
Faddeev-Popov operator 

Ar\x,y)Lo\y)^S,yJ2^AGt^\y)[^\y)-^''iy + e^)] 
+ Y.r'nA''^{y)uJ^{y + e,) 

C 

-Aliy~e,)u;^iy-e,)]}, (3) 
(x) = itr (K, u} [U^ix) + U^{x)+]) , 

with the generators of SU(2). As said above, we select 
always a random copy from the set of Gribov copies satis- 
fying the constraint in Eq. ([Ij through the minimization 
of the functional in Eq. thereby implementing the 
so-called minimal version of all of these gauges. 

The functional ([2]) can be minimized using the same 
methods as in Landau gauge [13] • In particular, a 
stochastic overrelaxation algorithm is used, analogous to 
the Landau-gauge case [2^. One should note that the 
required number of gauge-fixing sweeps increases with 
decreasing A. Hence simulations at smaller and smaller 
A become computationally more and more expensive. Al- 
though the number of gauge-fixing sweeps can be de- 
creased by multiple updates of each individual time-slice 
before changing the time-slice, this docs not decrease the 
total computing time needed, due to the increase of the 
number of updates on single time-slices. 

Thus, we limit our systematic investigations here to 
A > 10~^. Nevertheless, checks on small volumes down 
to A = 10"'^ do not show any qualitatively new effects. 

We distinguish furthermore between three types of 
Coulomb gauge fixing: 

• Simple Coulomb gauge, defined by setting A to 0. 

• Complete Coulomb gauge, defined by fixing the 
residual (perturbative) gauge freedom in the simple 
Coulomb gauge in a Landau-like manner Q. 

• Limiting Coulomb gauge, defined by the limit 
A — > 0. Note that while this limit is well-defined 
in the sense of a gauge condition, the Individ- 
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ual Green's functions do not necessarily have a 
pointwise-smooth liniit^. 

The last two versions of the Coulomb gauge should co- 
incide in the limit of A — )• @. We will determine the 
gluon propagator and the color- Coulomb potential for 
these three different prescriptions of fixing the Coulomb 
gauge. In all cases, the results coincide, as far as resolv- 
able. This will be discussed in more detail below. 

The details of the configurations employed for the de- 
termination of the propagator are given in Table H] Not 
all configurations have been used in the figures. The re- 
maining ones have been obtained to investigate various 
aspects, such as finite-size effects. In these cases, the re- 
sults are just stated in the text. The configurations used 
for the calculation of the color-Coulomb potential and 
associated quantities are hsted in Table HI 



III. PROPAGATORS IN A-GAUGES 

Setting A to a value different from 1 yields a gauge 
condition ([T|) that is no longer manifestly Lorentz (Eu- 
clidean) invariant. Thus, as in Coulomb gauge, the en- 
ergy po and the spatial momentum |p| have to be consid- 
ered as independent variables. 

Furthermore, the gluon is a vector particle. Hence, 
it is no longer possible to describe it by a single scalar 
function as in Landau gauge, but two independent ones 
are required. For practical purposes we choose these two 
independent functions to be^ 



-Doobo, IpI) 



1 



(N^ 1) 
1 



(4) 



< A^A'^ >, (5) 



where i and j run only over spatial indices. Z?''^ is eval- 
uated at zero d — I-momentum, and thus also at zero 
d-momentum, by setting the second term in the above 
projector to zero and multiplying by (d — 2)/((i — 1). 
This yields a unique and smooth limit from any direc- 
tion for vanishing d-momentum. Note that at A = 0, the 
second term in ([3]) is not contributing, since in this case 
PiAi = 0, due to the Coulomb gauge condition. The re- 
lation to the conventional choice in DSE studies [1^ is 
detailed in Appendix [X] 



^ However, our results suggest that the limit is smooth except at 
p = 0. 

^ Note that these are essentially the same definitions as were used 
in since there only the case po = has been investigated. 
This definition has been chosen for numerical convenience. See 
appendix lAlfor the relation to the definition of [521 ^^'^ ^ general 
decomposition. 



Herein, the gluon field in momentum space is ob- 
tained from the position-space variables by 



Alip) 



X 



where on a finite lattice the components P^^ of P have the 
integer values —N/2 -I- 1 , . . . , N/2 and the components 
Xfj, of X are the (integer) coordinates in the lattice rang- 
ing from to — 1 . The physical momenta p to which 
the integer lattice ones P correspond are 



Pi^ 



N 



The ghost propagator DcipoApl) Landau 
gauge, defined by the Fourier transform of the inverse 
Faddeev-Popov operator ([3]). As this operator is symmet- 
ric, the inversion can be done using a conjugate gradient 
method, as in ordinary Landau gauge (28j . To determine 
the renormalization effects in four dimensions, it is fur- 
thermore useful to define the dressing function G via 



Dg{poAp\) 



—^^—^G{po,\p\). 
Apo + P 



(6) 



Thereby the tree-level part l/{XpQ+p^) of the propagator 
is explicitly factored out. 



A. Propagators in 3 dimensions 

In three space-time dimensions, Yang-Mills theory in 
A-gauges is, as in Landau gauge, finite on the level of 
correlation functions. Hence no renormalization effects 
occur. In particular, the gauge-parameter A is not differ- 
ent from its tree-level value. 

It is worthwhile to study the three different propaga- 
tors, the one of the transverse gluon ([5]), the one of the 
temporal gluon (|3]), and the one of the ghost ([5]), indi- 
vidually. 

The transverse gluon propagator is shown for var- 
ious kinematic configurations and values of A in Figure [T] 
There are a number of interesting observations, outlined 
next. 

With respect to the spatial momentum, the clearly dis- 
tinct maximum of the transverse gluon propagator in 
Landau gauge moves towards larger momenta with de- 
creasing A, albeit rather slowly. As for the dependence 
on the temporal momentum, the dependence on po be- 
comes flat at low momenta and the flat region becomes 
more extended when decreasing A. However, the po — 
point is always lower, and may be influenced differently 
by finite-size effects. Already at A = 1/10, the maximum 
seems to be gone. The propagator at A = 1/100 is even 
more similar to the one of a massive particle. This is, 
however, a finite- volume effect. A comparison of different 
volumes, presented in Figure [21 shows that the volumes 
in Figure [1] are too small to resolve the position of the 
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FIG. 1: The transverse gluon propagator D*' for different values of A and different kinematic configurations, in the 3d case on a 
40^ lattice. The top-left panel is the propagator as a function of pure spatial momenta p, the top-right panel of pure temporal 
momenta po- The spatial momenta are measured along the a;-axis. Results shown are for A = 1, i.e. Landau gauge (circles), 
A = 1/2 (crosses), A = 1/10 (squares), A = 1/20 (triangles), A — 1/100 (stars) and A = (upside-down triangles). The lower 
panels show the full momentum dependence at A = 1 (middle- left panel), A = 1/10 (middle-right panel), A — 1/100 (bottom-left 
panel), and A = (bottom-right panel). All results are at /3 = 4.2. Also, A = corresponds to the simple- Coulomb-gauge 
case. Results at /? = 6.0 look similar, but are more strongly affected by finite-size effects. (Note that scales on vertical axes are 
different in individual sub-figures.) 



5 



TABLE I: Information on the configurations considered in our numerical simulations. The value "0 fixed" for A denotes the 
case when the residual gauge freedom (after setting A to 0) i s fixed according to 0|. The value of the lattice spacing a has been 
taken from [2^ for the three-dimensional case and from [sol. l3l| for the four-dimensional case. Sweeps indicates the number of 
sweeps [2^ between two consecutive gauge-fixed measurements. In all cases 200 thermalization sweeps have been performed. 
More details on the generation of configurations, error-determination, etc. can be found in [28l ]. 









d = 


3 dimensions 








d = 


4 dimensions 




A 


iV 


P 


a-' [GeV] 


Conf. 


Sweeps 


V^^" [fm] 


N 


13 


a-^ [GeV] 


Conf. 


Sweeps 


V^^" [fm] 


1 


40 


4.2 


1.136 


586 


50 


6.9 


22 


2.2 


0.938 


425 


50 


4.6 


1 


40 


6.0 


1.733 


689 


50 


4.5 


22 


2.5 


2.309 


532 


50 


1.9 


1 
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342 


50 


0.73 
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40 


2.2 


0.938 


50 




8.4 


1/2 


40 


4.2 


1.136 


466 


50 


6.9 


22 


2.2 


0.938 


380 


50 


4.6 


1/2 


40 


6.0 


1.733 


542 


50 


4.5 


22 


2.5 


2.309 


372 


50 


1.9 


1/2 














22 


2.8 


5.930 


362 


oO 


0.73 


1/10 


20 


4.2 


1.136 


481 


30 


3.5 














1/10 


20 


6.0 


1.733 


555 


30 


2.3 














1/10 


40 


4.2 


1.136 


424 


50 


6.9 


22 


2.2 


0.938 


311 


50 


4.6 


1/10 


40 


6.0 


1.733 


461 


50 


4.5 


22 


2.5 


2.309 


370 


50 


1.9 


1/10 














22 


2.8 


5.930 


399 


50 


0.73 


1/10 


60 


4.2 


1.136 


314 


70 


10 


40 


2.2 


0.938 


50 




8.4 


1/10 


60 


6.0 


1.733 


322 


70 


6.8 


70 


2.2 


0.938 


310 




14.7 


1/20 


40 


4.2 


1.136 


447 


50 


6.9 


22 


2.2 


0.938 


364 


50 


4.6 


1/20 


40 


6.0 


1.733 


463 


50 


4.5 


22 


2.5 


2.309 


307 


50 


1.9 


1/20 
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0.73 


1/20 
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50 
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1/20 














70 


2.2 


0.938 


336 




14.7 
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14 


2.2 


0.938 


418 


30 
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1/100 














14 
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420 


30 
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1/100 


40 


4.2 


1.136 


356 


50 


6.9 


22 


2.2 


0.938 


262 


50 


4.6 


1/100 


40 


6.0 


1.733 


478 


50 


4.5 


22 


2.5 


2.309 


367 


50 


1.9 


1/100 
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50 
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50 
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50 


4.5 


22 


2.5 


2.309 


340 


50 


1.9 

















22 


2.8 


5.930 


319 


50 
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fixed 
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377 


50 
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22 
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0.938 


434 


50 


4.6 


fixed 


40 


6.0 


1.733 


314 


50 


4.5 


22 


2.5 


2.309 


309 


50 


1.9 


fixed 














22 
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5.930 


401 


50 


0.73 



maximum at A ^ 1/10. It is visible that the maximum 
reappears at smaller momenta on a larger volume. Cor- 
responding calculations to track the maximum at even 
lower values of A become, however, exceedingly costly. 

For all three types of Coulomb gauge the transverse 
gluon propagators as functions of spatial momentum co- 
incide. This is not surprising. Quantities that are only 
defined on a given time-slice, such as the spatial gluon 
propagator, should be independent of how the inter-time- 
slice gauge degree of freedom is fixed. Therefore, in the 
figures only one case, that of the simple Coulomb gauge, 
is shown. 

We thus see that, in 3d Coulomb gauge, the transverse 
propagator for spatial momenta is similar to the one ob- 
tained in four dimensions in simple Coulomb gauge js^ . 
It has a maximum and is infrared-suppressed. 

On the other hand, the transverse propagator as a 



function of energy in Coulomb gauge is approximately 
constant. In the simple Coulomb gauge this is a con- 
sequence of the unfixed residual gauge degree of free- 
dom. That this persists even in the case of fully fixed 
Coulomb gauges is not trivial @, but has been deduced 
also in functional studies (ssj . In any case, in the alterna- 
tive ways to fix the remaining gauge freedom, the gauge 
condition cannot impose a strict constraint on the gluon 
fields on any given configuration'^. Hence, the fields still 
wash out when averaged over configurations, leaving only 
an instantaneous propagator. For this reason, the prop- 



^ Consider, e.g., 9o^O = 0. Such a constraint is only ful- 
filled on average, < doAo >= 0. The integrated constraint 
do J d'^~^XiAo{xo,Xi) = is, however, true for any individual 
configuration. 
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TABLE II: Configurations used for the determination of the 
color- Coulomb potential and the residual-gauge-symmetry or- 
der parameter. For details, see Table (T) All calculations have 
been performed in four dimensions. The columns with various 
/3's indicate the number of configurations used. The number 



of sweeps 


was 


always 50. 


















V ' [fm] 


A 


N 


/3 = 2.2 


13 = 2.5 


(3 = 2.8 


at ^ = 2.2 


1 


4 


188 


188 


188 


0.84 


1 


6 


190 


190 


190 


1.3 


1 


8 


191 


191 


191 


1.7 


1 


10 


192 


192 


192 


2.1 


1 


12 


193 


193 


193 


2.5 


1 


16 


195 


195 


195 


3.4 


1 


20 


196 


196 


196 


4.2 


1 


24 


104 




197 


5.0 


1/10 


4 


188 


188 


188 


0.84 


1/10 


6 


190 


190 


190 


1.3 


1/10 


8 


191 


191 


191 


1.7 


1/10 


10 


192 


192 


192 


2.1 


1/10 


12 


193 


193 


193 


2.5 


1/10 


16 


195 


195 


195 


3.4 


1/10 


20 


180 


196 


196 


4.2 


1/10 


24 


94 


298 


197 


5.0 


1/100 


4 


188 


188 


188 


0.84 


1/100 


6 


190 


190 


190 


1.3 


1/100 


8 


191 


191 


191 


1.7 


1/100 


10 


192 


192 


192 


2.1 


1/100 


12 


193 


193 


193 


2.5 


1/100 


16 


100 


195 


195 


3.4 


1/100 


20 


40 


192 


138 


4.2 


1/100 


24 


29 


197 


197 


5.0 





4 


188 


188 


188 


0.84 





6 


190 


190 


190 


1.3 





8 


191 


191 


191 


1.7 





10 


192 


192 


192 


2.1 





12 


193 


193 


193 


2.5 





16 


195 


195 


195 


3.4 





20 


115 


196 


196 


4.2 





24 


70 


197 


197 


5.0 


fixed 


4 


188 


188 


188 


0.84 


fixed 


6 


190 


190 


190 


1.3 


fixed 


8 


191 


191 


191 


1.7 


fixed 


10 


192 


192 


192 


2.1 


fixed 


12 


193 


193 


193 


2.5 


fixed 


16 


195 


195 


195 


3.4 


fixed 


20 


82 


196 


196 


4.2 


fixed 


24 


90 


197 


197 


5.0 



agator is constant as a function of the energy. This is 
observed in all three types of Coulomb gauge and shown 
in Figure [T] for the simple-Coulomb case. Note that this 
argument does not apply to the zero-energy mode. This 
mode coincides vifith the zero mode depending on spatial 
momenta, leading to the same value at zero d-momentum 
for the propagator, irrespective of the dependence on spa- 
tial momentum or on energy. If the propagator as a func- 
tion of spatial momentum thus vanishes in the infinite- 
volume limit, it may have a discontinuity at zero as a 



Volume dependence of D"^ as a function of p 




FIG. 2: The transverse gluon propagator D*"^ as a function of 
po at A = 1/10 for different volumes, in the 3c? case. Open 
and full symbols correspond to j3 = 6.0 and /3 = 4.2, respec- 
tively. The lattice sizes are 20^ (circles), 40'^ (squares), and 
60^ (triangles). 



function of energy. 

Concluding, the spatial propagator D*^'^{p) is instan- 
taneous and infrared-suppressed in all versions of the 
Coulomb gauge considered here. 

Although not shown, a further observation is that for 
the spatial propagator as a function of p the dependence 
on /3 at fixed lattice size weakens with decreasing A. How- 
ever, at the same time the /3-dependence increases for 
momenta along the time direction. This can indicate 
either a change in the volume dependence of the propa- 
gators with a change in A, or a different /3-dependence. 
The former is more likely, as can be seen in Figured The 
volume dependence is stronger than in correspon ding cal- 
culations in Landau gauge in three dimensions |28l. [29j. 
It is furthermore visible that the ^-dependence dimin- 
ishes with increasing volume, making it more probable 
that this is a finite- volume effect, rather than a finite-/? 
effect. 

The case of the temporal propagator is at the first 
glance a bit simpler, as Dqq{pq, 0) vanishes identically for 
Po 7^ due to the gauge condition. The remaining results 
are shown in Figure [S As a function of spatial momenta, 
the maximum moves towards smaller momenta with de- 
creasing A. Again, already at A = 1/10, resolving the 
maximum is no longer possible for these volumes. At 
the smallest values of A, the propagator looks more like 
a divergent function, rather than an infrared-suppressed 
function, except at the zero-momentum point for finite 
A. However, this is again only a finite- volume effect, 
as is visible in Figure ID This changes drastically for 
the Coulomb gauges, where the propagator diverges for 
p 0, independently of po, as it does already in pertur- 
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FIG. 3: The temporal gluon propagator Dqq for different values 
of A and different kinematic configurations, in the 3d case on a 
40"^ lattice. The top-left panel is the propagator as a function of 
pure spatial momenta p. The momenta are measured along the 
X-axis. Results shown are for A = 1, i.e. Landau gauge (circles), 
A = 1/2 (crosses), A = 1/10 (squares), A = 1/20 (triangles), A = 
1/100 (stars), and A = (upside-down triangles). The panels 
below show the full momentum dependence at A = 1 (middle- 
left panel), A = 1/10 (middle-right panel), A — 1/100 (bottom- 
left panel), and A = (bottom-right panel). All results are at 
(3 — 4.2. Also, A = corresponds to the simple-Coulomb-gauge 
case. Results at /3 = 6.0 look similar, but are more strongly 
affected by finite-size effects. (Note that scales on vertical axes 
are different in individual sub-figures.) 
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bation theory. Thus, for all three definitions of Coulomb 
gauge the infrared-divergent gluon propagator is recov- 
ered. Note in particular that also the value at p = is 
non-vanishing in the Coulomb limit. Furthermore, the 
limit A — > is thus discontinuous for vanishing spatial 
momentum. It is, however, smooth for any non-vanishing 



spatial momentum. Concluding, the propagator is found 
to be instantaneous in all three Coulomb-gauge defini- 
tions. 

Finally, the ghost propagator is shown in Figure [S] 
The dependence on A is clearly seen to increase with 
a decreasing ratio of po/l??!- There is essentially no A- 
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FIG. 4: The temporal gluon propagator Doo as a function of 
IpI at A = 1/10 for different volumes, in the 3d case. Open 
and full symbols correspond to /? = 6.0 and /3 = 4.2, respec- 
tively. The lattice sizes are 20'^ (circles), 40'^ (squares), and 
60^ (triangles). 



dependence at po = 0. Furthermore, as in Landau gauge, 
always a strong infrared enhancement is visible. Hence 
the most interesting case is the one of pure temporal mo- 
menta. In this case the propagator seems (to leading 
order) to be scaled by a constant factor, which is both 
A- and /3-dcpcndcnt. This is tested explicitly in Figure 
[SI Multiplying the dressing function by A", with a cho- 
sen to be 0.15 for 13 = 4.2 and 0.1 for = 6.0, yields 
roughly agreeing dressing functions. This is a genuine /S- 
depcndent effect, as the same exponents appear on a 20'^ 
lattice at the same values of /3. This may indicate that 
for this correlation function the scaling regime is moved 
as a function of /3. The reason for this movement of the 
scaling region could be that the theory in Coulomb gauge 
is no longer finite, and renormalization has to be taken 
into account, even in three dimensions. If this should oc- 
cur smoothly, the scaling region without renormalization 
should move to larger (3 with decreasing A. An alternative 
explanation of this may be a combination of finite-size ef- 
fects, leading seemingly to a /3-dependencc, and a genuine 
non-trivial dependence of the dressing function G on A. 
Such a non-trivial dependence is in fact found in func- 
tional calculations in the far infrared (see Appendix \K\ 
and [2^). However, if this is the correct interpretation, 
than this non-trivial A-dcpendencc seems to pertain to 
rather large momenta. 

If the exponent a should be non-zero in the infinite- 
volume and continuum limit, then this would imply that 
the ghost propagator diverges for all po, independently of 
the spatial momentum, in the limit A — >■ 0. This would be 
the behavior expected in the Gribov-Zwanziger scenario 

am. 



It is not possible with the method used here to de- 
termine the ghost propagator in the case of A = 0, as 
the corresponding Faddeev-Popov operator is different 
from ([3]). In this case, it would be necessary to invert 
the Faddcev-Popov operator in the sub-space orthogonal 
not onlyto constant gauge transformations, as in Landau 
gauge [28| , but also to all purely temporal gauge transfor- 
mations. If, however, the limiting Coulomb gauge A — 0, 
the simple and the complete Coulomb gauges with A = 
coincide also for the ghost propagator (as in the case 
of the gluon propagator), then an infrared-divergent, in- 
stantaneous ghost propagator would be obtained. This 
is the limiting behavior that is suggested from Figure [Sj 

Concluding, the results show how the propagators 
evolve smoothly with A from a Landau-like behavior to 
a Coulomb-like behavior. However, below a certain mo- 
mentum, which decreases with decreasing A, the prop- 
agators always show qualitatively a Landau-like behav- 
ior. In particular, the temporal propagator as a func- 
tion of spatial momenta shows a more and more diverg- 
ing behavior, as in Coulomb gauge, at intermediate mo- 
menta, before becoming infrared-suppressed, as in Lan- 
dau gauge. On the other hand, the transverse propaga- 
tor as a function of transverse momenta becomes more 
and more infrared-suppressed, as it is in Coulomb gauge. 
Finally, the ghost propagator as a function of pure tem- 
poral momenta shows a rather pronounced dependence 
on A and /3. 

All in all, the propagators show a smooth development 
towards a Coulomb-like behavior in an intermediate mo- 
mentum regime. Only in the far infrared does the original 
Landau-like behavior persist. The limit A ^ thus cor- 
responds to the limit in which the onset momentum of 
Landau-like behavior is moved towards zero momentum. 
As the Coulomb-like behavior and the Landau-like be- 
havior are different, this limit is not smooth on the level 
of correlation functions. These results agree qualitatively 
with the naive momentum- rescaling argument of @. 

An interesting observation is the coincidence of the 
three different definitions of the Coulomb gauge. While 
the simple Coulomb gauge leaves the component ^0 com- 
pletely unfixed, the limiting and the complete Coulomb 
gauge fix this component. However, in both cases the 
gluon propagators coincide within statistical errors. That 
is a consequence of the additional constraint not being a 
local condition on the fields. It only ensures OqAq = 
on average, just as in the case of an unfixed temporal di- 
rection. To completely understand the properties of the 
limit of the interpolating gauge, this has to be investi- 
gated further, but is already in line with arguments from 
functional calculations 1331. 



B. Propagators in 4 dimensions 

The situation in four space-time dimensions turns out 
to be qualitatively the same as in three dimensions. Still, 
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FIG. 5; The ghost dressing function G for different values of A 
and different kinematic configurations, in the 3d case on a 40'^ 
lattice. The top-left panel is the propagator as a function of 
pure spatial momenta p, the top-right panel of pure temporal 
momenta po- The spatial momenta are measured along the x- 
axis. Results shown are for A = 1, i.e. Landau gauge (circles), 
A = 1/2 (crosses), A = 1/10 (squares), A = 1/20 (triangles), and 
A — 1/100 (stars). The lower panels show the full momentum 
dependence at A = 1 (middle- left panel), A = 1/10 (middle-right 
panel), and A = 1/100 (bottom- left panel). All results are at 
P = 4.2. Also, A = corresponds to the simple-Coulomb-gauge 
case. Results at /? = 6.0 are essentially identical. 



quantitatively, differences appear. One is due to the 
fact that finite- volume effects are harder to control, since 
calculations arc in general more costly. The second is 
the need for renormalization even outside the Coulomb 
gauge: The theory is no longer finite in four dimensions. 
This point makes the situation a bit more intricate, as the 
gauge parameter A also gets renormalized by a renormal- 
ization constant Z\ . This influences of course the defini- 
tion of the ghost dressing function (jH]). In addition, the 
ghost dressing function itself gets renormalized by the 
wave- function renormalization constant Z^. The corre- 



sponding relative renormalization constants between two 
/3-values will be determined by 

ZsiPi) ^ gG(o,b1,/?i) . . 

Z^{(3,) Z?G(0,b1,/32) ^' 

ZxiM DG{po,0,P2)DGiOAp\,M . . 

ZAfi2) ^ DG(po,0,;3l)DG(0,b1,/32)' ^' 

The absolute normalization can be chosen at will, and is 
of no importance here. The second relation is based on 
(pl , and would be exact only if the dressing function did 



10 



Ghost dressing function G as a function of p 



> 

0) 

O 



B 3 



# 

o 

■ 

□ 

■ 

A . 



-B 
A 



0.5 



1.5 



Pp [GeV] 



FIG. 6: The rescaled ghost dressing function X"G for different 
values of A for temporal momenta po, in the 3d case. Symbols 
used are the same as in the top row of Figure O with open 
symbols corresponding to P — 6.0 and closed to /? = 4.2. 



not depend on /?. Within the domain of perturbation 
theory, this dependence exists, but is only logarithmic. 
For the quite similar values of (3 employed here, relation 
([5]) is hence an acceptable approximation at sufficiently 
large momenta. 



TABLE III; Ratios ([THH]) of the renormalization constants. 
The matching condition was a continuous matching at the 
largest edge-momenta of the smaller /3-value. Note that at 
A = 1 the gauge parameter is not renormalized because Lan- 
dau gauge is a fixed-point of the renormalization group. At 
A = no results on the ghost propagator are available, see 
text. At A = 1 the two gluon renormalization constants coin- 
cide trivially. Note that, at least in the simple Coulomb gauge, 
it is strictly speaking not possible to renormalize purely mul- 
tiplicatively [3g |. 



A 


Pi 


/?2 


Z3(M 






Z3(/32) 






1 


2.2 


2.5 


1.12 


1 


1.23 


1 


2.5 


2.8 


1.08 


1 


1.28 


0.5 


2.2 


2.5 


1.12 


0.95 


1.23 


0.5 


2.5 


2.8 


1.09 


0.95 


1.28 


0.1 


2.2 


2.5 


1.12 


0.65 


1.23 


0.1 


2.5 


2.8 


1.10 


0.91 


1.28 


0.05 


2.2 


2.5 


1.12 


0.55 


1.23 


0.05 


2.5 


2.8 


1.11 


0.83 


1.28 


0.01 


2.2 


2.5 


1.12 


0.4 


1.23 


0.01 


2.5 


2.8 


1.12 


0.77 


1.28 





2.2 


2.5 






0.65 





2.5 


2.8 






0.74 



In case of the gluon propagator, the corresponding ra- 



tios of renormalization constants will be defined by 

. i?'^(0,b1,/3i) 



zim 



^*'-(o,Ip1,/32) ' 

^00(0, b1,/?i) 
^00 (0, b1,/32) 



(9) 
(10) 



Note that due to the lack of Euclidean invariance both 
propagators have to be renormalized separately [2], 
and therefore independent wave- function renormalization 
constants Z3'' for the transverse gluon propagator and 
for the temporal gluon propagator have to be introduced. 
This is in contrast to Landau gauge. Also, note that 
when reaching Coulomb gauge the combination q^Dqq is 
a renormalization group invariant, and the corresponding 
inverse renormalization constant of the propagator alone 
yields the renormalization of the coupling constant 

The approximate renormalization constants are given 
in Table Hill The ratios for the transverse and the tem- 
poral part of the gluon propagator are found to be essen- 
tially equal, and thus only one is listed. Of course, these 
ratios cannot be determined precisely, as the same mo- 
menta are not available for different lattices. The renor- 
malization has been performed to obtain a smooth con- 
nection at the largest edge-momentum of the lattice with 
the larger physical volume at the same lattice volume^. 
In addition, statistical fluctuations and systematic errors 
also prevent a precise measurement within the present 
scope, so these results arc more qualitative than quan- 
titative. However, the error on the ratio is of the same 
size as the statistical and systematic errors of the propa- 
gators, and thus a few, up to ten, percent for statistical 
errors and similarly for the systematic errors. The abso- 
lute normalization is chosen not to change the results at 
P = 2.5. 

These results indicate that the wave-function renor- 
malization constants do not depend markedly on A. How- 
ever, the renormalization constant of A itself depends 
on A. This was somewhat anticipated. When moving 
further away from a fixed point of the renormalization 
group, renormalization effects should become stronger 
and stronger with "distance" . That Landau gauge is a 
fixed point is explicitly demonstrated by the decrease in 
Z\ with increasing /3. From here on, only renormalized 
results will be shown. 

Aside from these renormalization effects, the results 
in four dimensions show the same qualitative behavior 
as in three dimensions. However, this is only true when 
comparing roughly the same volumes, or, in fact, slightly 
smaller volumes in the three-dimensional case, to the vol- 
umes available in four dimensions. In particular, nearly 
none of the infrared effects visible in three dimensions can 



^ The ultraviolet properties of the gluon propagator are only 
weakly affected by finite-volume effects. Hence the difference of 
physical volume can be ignored for the purpose of determining 
the renormalization constants at ultraviolet momenta. 
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(yet?) be seen at the volumes used in four dimensions. 
One important exception below will be the transverse 
gluon propagator as a function of transverse momentum 
at very small A. This has already been discussed previ- 
ously in 0. Hence, here only a more compact presenta- 
tion of the results will be given compared to the one in 
three dimensions in the previous section. In particular, 
as the results for non- vanishing energy and non- vanishing 
spatial momenta only interpolate smoothly between dif- 
ferent behaviors if one of them is zero, here only these 
cases will be shown. Of course, in these kinematic con- 
figurations the results also look qualitatively similar to 
the three-dimensional case. 

The results for the transverse gluon propagator are 
shown in Figure [T] As a function of spatial momen- 
tum \p\ , the propagator becomes more and more infrared- 
suppressed, as already observed for three dimensions. No 
clear maximum is visible at large A, owing to the small 
volumes in four dimensions. However, for the large vol- 
umes at values of A < 0.1, a shallow maximum is ob- 
served. This maximum moves to larger momenta with 
decreasing A, as in the three-dimensional case. This can 
also be seen in a direct comparison of different volumes, 
shown in Figure [51 This confirms the results obtained 
in 0. The observed infrared-suppressed gluon propaga- 
tor in Coulomb gauge is in accordance with the results 
of 0, [13, HI] , though here the momentum-dependent as- 
pects of renormalization in Coulomb gauge still have to 
be taken into account (36j . 

The behavior as a function of temporal momentum po 
is also reminiscent of the results in three dimensions. But 
due to the restriction in volume the effects are even less 
pronounced. Still, it is visible how the propagator be- 
comes flatter with decreasing A and is constant at A = 0, 
except at zero momentum. 

Hence, up to the effects of the smaller volumes and 
renormalization, the transverse gluon propagator be- 
haves qualitatively as in three dimensions. In particu- 
lar, as a function of spatial momentum, it becomes more 
infrared-suppressed with decreasing A. 

The results for the temporal gluon propagator are 
shown in Figure [5] The situation is essentially the same 
as in three dimensions. The propagator as a function of 
spatial momenta becomes more and more enhanced in 
the infrared. Of course, for the volumes accessible here, 
it is not possible to check whether it then bends over 
again in the far infrared, as it was observed in three di- 
mensions in sufliciently large volumes. When adding a 
dependence also on the temporal momenta, it is found 
that the propagator increases with decreasing A as well. 
Note that this implies a much stronger falloff towards 
zero temporal momenta, which always occurs due to the 
gauge condition for vanishing spatial momenta. 

The dependence on volume is less pronounced (and 
less spectacular so far) than in case of the spatial gluon 
propagator, and hence will not be shown explicitly. 

Finally, the ghost dressing function is shown in Fig- 
ure 1101 It shows a clearly infrared enhanced behavior at 
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FIG. 7: The transverse gluon propagator D^'^ for different 
values of A. The top panel is the propagator as a function of 
pure temporal momenta po, the bottom panel of pure spatial 
momenta p. The spatial momenta are measured along the x- 
axis. Results shown are for A = 1, i.e. Landau gauge (circles), 
A = 1/2 (crosses), A = 1/10 (squares), A = 1/20 (triangles), 
A = 1/100 (stars), and A = (upside-down triangles). All 
results are at /3 = 2.2. Also, A = corresponds to the simple- 
Coulomb-gauge case. Results at larger /? look similar, but 
are more strongly affected by finite-size efi'ects. In case of the 
spatial momenta, the results for A = 1 come from a lattice of 
size 40", while A = 1/10, A = 1/20, and A = 1/100 are from a 
lattice of size 70'*. In all other cases, the lattice size was 22*. 



small momenta, essentially unaltered compared to Lan- 
dau gauge. Only in the pure temporal case is an addi- 
tional increase with A for all momenta observed which is 
not a pure renormalization effect, but instead similar to 
the situation in three dimensions. 

As in three dimensions, it is possible to let all curves 
at p = collapse by multiplying with A°. In contrast 
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Spatial gluon propagator d"^ as a function of |^| and volume 
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FIG. 8: The transverse gluon propagator D'"^ as a function 
of IpI at A = 1/100 for different volumes, in the 4d case. 
Circles, squares, triangles, stars, and open circles correspond 
to systems of size 14*, 22*, 32* (from 0), 40*, and 70*, all at 
/3 = 2.2. 



Temporal gluon propagator D as a function of |p| 
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FIG. 9: The temporal gluon propagator Doo for different val- 
ues of A as a function of pure spatial momenta p, in the 4d 
case. The momentum is measured along the x-axis. Symbols 
have the same meaning as in Figure [T] 



to three dimensions, the exponent a is essentially /3- 
independent, as a consequence of renormalization. It in- 
creases slowly with A, more or less leveling off at the 
smallest achieved A value at around a « 0.08. This 
rescaling leads to an apparent weaker divergence of the 
ghost dressing function. However, with such a small a, 
the effect is rather weak. 

Thus the results in four dimensions arc, up to finite- 
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FIG. 10: The ghost dressing function G for different values of 
A as a function of pure spatial momenta p, in the 4d case, in 
the top panel. The dependence on pure temporal momenta po 
is given in the bottom panel. Results shown are for A = 1, i.e. 
Landau gauge (circles), A = 1/2 (crosses), A = 1/10 (squares), 
A = 1/20 (triangles), and A = 1/100 (stars), all at /3 = 2.2 on 
a 22* lattice. Results at larger /? are essentially identical. 



volume effects, qualitatively the same as in three dimen- 
sions. When decreasing A then at least at mid-momenta 
a Coulomb-like behavior is observed. In particular the 
transverse gluon propagator as a function of spatial mo- 
mentum becomes more infrared-suppressed, the tempo- 
ral propagator as function of spatial momentum becomes 
more mid-momentum enhanced. The ghost seems to be 
more or less inert, but as it is the quantity usually least 
affected by changes which are effectively volume changes, 
this is not surprising. 

For Coulomb gauge the same applies, as was already 
said in the three-dimensional case: On the level of the 
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FIG. 11: The ratio of the transverse gluon propagator ([5]) to 
the instantaneous would-be Coulomb-gauge propagator (lllf) 
as a function of spatial momentum for A = 1 (circles), A — 
1/10 (squares), A — 1/20 (triangles), and A = 1/100 (stars). 
The lattice was of size 40'' at 13 = 2.2. 



gluon propagators all definitions of the Coulomb gauge 
are indistinguishable. This approaeh can be made even 
quantitative, e. g. by comparing the gluon propagator as 
a function of spatial momentum defined according to ([5]) 
with the would-be Coulomb-gauge propagator, i.e. the 
time-average of (O 
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PzPj 
^2 
P 



<Ant,p)A-{t,-p) > 

{N^-l)id~2)V 



(11) 

This comparison is made in Figure 1111 It is explicitly 
visible how the transverse gluon propagator tends to the 
instantaneous one, in particular in the infrared. In fact, 
at the smallest non-vanishing momenta the rate of ap- 
proach is roughly proportional to the unrenormalized A 
beyond Landau gauge. 



IV. COULOMB POTENTIAL 

One remarkable property of Coulomb gauge is the ex- 
istence of an upper bound to the conventional Wilson 
potential in terms of a gauge-dependent quantity, the 
so-called color-Coulomb potential This condition 

is often called informally the "no-confinement-without- 
Coulomb-confinement" one, and has been repeatedly 
studied [i,[li,|23. 

The Coulomb potential can be straightforwardly de- 
fined by use of Unk variables as 



VciR,0)^Eo 



In 



--tr{Uo{x,t)U+{y,t)) 



" — y\ 

(12) 



FIG. 12: The color-Coulomb potential (|12p as a function of 
distance on a 24* lattice at /3 = 2.2. Other volumes and 
discretizations follow the expected pattern. The gauges are 
A = 1 (circles), A = 1/10 (squares), A = 1/100 (stars), sim- 
ple Coulomb gauge (upside-down-triangles) and fully fixed 
Coulomb gauge (triangles). 



where Eq is a constant energy offset, which can be re- 
moved by renormalization of the potential. For conve- 
nience, it is here always chosen such that the Coulomb 
potential is zero at zero distance. Of course, this identi- 
fication is only possible up to ©(a^), but studies at vari- 
ous /? indicated that the effects are small, but discernible. 
However, they affected essentially only the value of the 
string tension, but neither the functional shape of Vc, 
nor its dependence on A, and an extensive analysis is 
therefore skipped. 

The result is shown in Figure [T^ It is clearly vis- 
ible that the color- Coulomb potential for all gauges, 
except Coulomb gauge, is asymptotically flat. How- 
ever, with decreasing A, a linear rise is seen initially at 
small distances. Finally, all three possible definitions of 
the Coulomb gauge, simple Coulomb gauge, fully-fixed 
Coulomb gauge, and the limit A coincide in case of 
the color-Coulomb potential. 

This result corroborates therefore the finding that cor- 
relation functions in interpolating gauges at small dis- 
tances behave Coulomb-gauge-like, and only at large dis- 
tances revert back to a Landau-gauge behavior. Hence, 
at short range, the interpolating gauges are almost 
Coulomb-like. The notion of small is therefore neces- 
sarily set to be of order 1/(AAym), though, of course, 
the functional shape could be different. This is only an 
estimate, but in line with the interpretation given in 0. 

The fact that once more all three implementations of 
Coulomb gauge are the same is at this point not surpris- 
ing. However, there exists an order parameter Q for the 
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FIG. 13: The remnant-residual-symmetry order parameter Q 
(|13|l as a function of L = aNx at (5 = 2.2. Though the pre- 
factors appearing are strongly /3-dependent, as investigations 
ai P — 2.5 and /3 — 2.8 have sliown, the functional behavior 
remains the same. The gauges are A = 1 (circles), A = 1/10 
(squares), A = 1/100 (stars), simple Coulomb gauge (upside- 
down-triangles) and fully fixed Coulomb gauge (triangles). 



residual gauge freedom of Coulomb gauge, defined as Q 



Q 

u{t) 




(13) 



N.,NyN, 



If for N^NyNzNt — > oo Q is vanishing, the residual gauge 
symmetry is unbroken. Since it has been shown that this 
is equivalent to a non-zero Coulomb string-tension Q , it 
is possible to anticipate from Figure [T^] already the re- 
sult that Q will vanish for all definitions of the Coulomb 
gauge. In fact, this is seen in the results presented in 
Figure [131 In all definitions of the Coulomb gauge Q 
vanishes with volume. This also implies that although 
the residual gauge symmetry has been fixed in the fully- 
fixed version of Coulomb gauge and by taking the limit of 
the interpolating gauge, the residual gauge symmetry is 
not broken on average. On the other hand, for all gauges 
except Coulomb gauge it is once more visible how for de- 
creasing A at small volumes first a Coulomb-like pattern 
emerges before the original Landau-gauge pattern, with 
broken residual gauge symmetry, becomes manifest. 

It should be noted that^ Q = implies U{t) = 0, i. e., 
a zero matrix. The second type of Coulomb gauge, in 
which the temporal gauge freedom of Coulomb gauge is 



^ We are grateful to Daniel Zwanziger for pointing this out. 



fixed in a Landau-like manner, is actually the condition 
to minimize 

tr t/o(a;,<) cx ^tr C7(i). 

x^y^z,t t 

In the confining case, this sum is thus always exactly zero 
in the infinite- volume limit, and therefore the correspond- 
ing gauge-fixing step is trivial. Thus, in the infinite- 
volume limit the fully-fixed Coulomb gauge and the sim- 
ple Coulomb gauge coincide. That this appears to be 
also the case at finite volume, as found here, appears 
non-trivial. 



V. SUMMARY 

Summarizing, a systematic investigation of correlation 
functions of Yang-Mills theory in the class of A-gauges in- 
terpolating between the Coulomb and the Landau gauge 
in three and four dimensions has been performed. 

The main result is threefold. On the one hand, 
the gauge parameter interpolating between both gauges 
parametrizes a certain length scale. For distances shorter 
than this scale, all the determined correlation functions 
show an essentially Coulomb-like behavior, at least up 
to momenta where perturbation theory becomes the 
dominant contribution. For distances longer than this 
scale, the behavior is essentially that as in the Landau 
gauge. There is no indication that this separation of 
Coulomb-like and Landau-like behavior will change for 
larger volumes and/or finer discretizations, even though 
the asymptotic infrared behavior in Landau gauge is still 
sensitive to such changes. 

Furthermore, the investigation of different realizations 
of the Coulomb gauge, with and without fixing the resid- 
ual gauge degree of freedom, did not show any effect for 
any of the correlation functions. This is in accordance 
with the reasoning of @ and of the continuum investiga- 
tions in j33| . 

Finally, the evolution of correlation functions with this 
gauge parameter from Landau gauge to Coulomb gauge 
is smooth, with the only possible exception of the point 
at zero four-momentum for the temporal gluon propa- 
gator: In the interpolating gauge it is forced to vanish 
there, while in Coulomb gauge it is expected to diverge. 
However, for any non-zero four-momenta, its behavior is 
smooth. 

Thus the interpolating gauge indeed turns into an in- 
terpolation of scales. The results interpolate in the far 
infrared from a Landau-like behavior over a Coulomb-like 
behavior to the pcrturbative A-likc behavior. Hence the 
name of interpolating gauges is truly justified. 

Concluding, it is remarkable that all energy- 
dependence is actually vanishing in all definitions of the 
Coulomb gauge. Thus, all investigated correlation func- 
tions become time- independent, and therefore do not de- 
scribe propagating degrees of freedom, in contrast to 
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physical (bound) states. In this sense, Coulomb gauge 
can be regarded as a physical gauge, for any of the inves- 
tigated versions. 
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Appendix A: Infrared analysis using 
Dyson-Schwinger equations 



In this appendix, connection to the commonly used 
framework for interpolating gauges in functional meth- 
ods will be made. Furthermore, it will be shown that 
the persistence of the Landau-like behavior in the far in- 
frared (at least for non-pcrturbativc A-gaugcs) allows for 
a scaling solution not only in four dimensions, as shown 
in (22| . but also in three dimensions. 

In functional calculations, it is convenient to use in- 
stead of the quantities Dqq and D*'', defined in equations 
(|4]-[5]) above, different ones. Writing the gluon propagator 
as (no summation implied) 



doo and d*'' are related to the ones defined in (|4H5]) by 



(Al) 



defines two scalar functions, doo and d*"^. 

Due to the appearance of \p\ in the denominator, this 
possibility to parametrize the gluon propagator is not 
too useful in lattice calculations. The scalar functions 



D, 



00 



pI 



p^{pI+p^) 

Pt 

p^{pI+p^) 



Dn 



The main difference in the results between both defini- 
tions is that the maximum in spatial (temporal) direc- 
tions^ is more pronounced for d''' (-D''^) than for Z?''' (d*''), 
when comparing both results with the lattice data. For 
the results presented in the main text, this is only a weak, 
quantitative effect. 

The DSEs and their solutions can then be obtained 
along the same lines as in four dimensions [22|. The 
equations are given by 

DcipT^ = PP (A2) 
- J d^kPt^Df,,,{k)p^DGip + k) 



D^,{p) = D„{p)D%{p)D„,{p) 



d^kD^p{p)kpkcrD^^{p) 



(27r)3 

X DG{k)DG{p + k), 
PfL = Hpp = {po\,p)^ , 

using the metric H = diag(A, 1, 1), and tl denotes a tree- 
level quantity. This is already an approximate system us- 
ing the same truncation as in [2^ : Only tree-level terms 
and the self-consistent leading terms for a scaling solu- 
tion, i.e., terms with at least one ghost-line, are retained. 
Furthermore, a pcrturbativc color structure is assumed. 

As in four dimensions [l^, equations (jA2HA3[) can be 
transformed so that they have the same form as in Lan- 
dau gauge. Since the gauge parameter is not renormal- 
izcd in three dimensions, it is sufficient to introduce the 
quantities 



Xd^p 
Dp,{p) 

DGiP) 



d^p 

Dg{p)- 



Note that this implicitly defines in general rather compli- 
cated functions D^^i, and Dg- Rewriting equations (jA3l - 



® Spatial directions implies still a non-zero temporal component, 
as otherwise both definitions coincide. 
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H) in terms of these new variables leads to 



Dcip) ^ = PP 



(A3) 



(27r)3 
Ca9^ 



D^p{p)D%{p)D,,[p) 



(A4) 



d? kDfj,p (p) kp ka Day (p) 



(2^)3 J - -^f^P 
X DG{k)DG{p+~k). 

Except for the tree-level terms, this system of equations 
looks now formally as in Landau gauge. In the gluon 
equation, (jA4[) . the tree- level term will turn out to be 
sub-leading, and can thus be neglected. In case of the 
ghost equation, this is in general more subtle. However, 
when making an infrared ansatz of type 



DaiP) 



P 



2k-2 



Sp.-^]Dpyip) ^ f'-^ 



the system needs to be renormalized to possess a solution. 
It is then possible to argue as in four dimensions [l^l that 
the tree-level term can be removed in the renormalization 
process. Thus the system becomes completely equivalent 
to the Landau gauge case, and thus has the same two 
solutions ySl 



2k + i = - 



1 

2 
1 

^2 

-0.39. 



(A5) 



In particular, this leads to the infrared behavior for the 
original functions of 



^\p) 
doo{p) 



1 




+ A2p2)l-K 




A 








AV 




+ A2p2)2-t 



(A6) 



The ghost diverges at zero p^ . Its anomalous scaling with 
A for p = is likely directly related to the appearance 
of the factor A^^^k jj^ (|A6p . The transverse gluon propa- 
gator is consistently suppressed from all directions. The 
temporal gluon propagator vanishes identically for p = 0. 
From all other directions it is also infrared-suppressed, 
but this is seemingly not so much due to the exponent t, 
but due to the explicit factor of "p^ in the numerator. In 
fact, 2 — t is not larger than zero. 
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